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TECHNICAL NOTE D-859 

UNSTEADY AERODYNAMIC FORCES ON A SLENDER BODY OF 
REVOLUTION IN SUPERSONIC FLOW 
By Reuben Bond and Barbara B. Packard 


SUMMARY 


Linearized slender-body theory is applied to the computation of 
aerodynamic forces on an oscillating, or deforming, body in supersonic 
flow. The undeformed body is a body of revolution and the deformed body 
is represented by movement of a line through the centers of the cross 
sections which are assumed to remain circular. The time dependence is 
based on sinusoidal motion. 

For a body of vanishing thickness the slender-body theory yields 
the apparent mass approximation as it is obtained for incompressible 
crossflow around a cylinder. 

Both linearized slender-body theory and the apparent mass approxi- 
mation are used to calculate the pitching-moment coefficients on a rigid 
slender body with a parabolic arc nose cone, and these coefficients are 
compared with some experimental results. 


INTRODUCTION 


The problem of predicting the dynamics of launch vehicles requires 
knowledge of the aerodynamic forces on an oscillating and/or a deforming 
body in high-speed flight. These forces are needed as an input to both 
the dynamic stability problem and as an input into the aero-servo-elastic 
problem. Especially important is the part of the ascent where the velocity 
is high and the atmospheric density is still great enough to cause signif- 
icant aerodynamic reactions. The present analysis is concerned with this 
region which is generally the low supersonic part of the flight. This 
problem is considered in chapter 12 of reference 1 where references are 
given to other work. None of the references consulted have shown the 
type of slender-body approximation used here. 

The analysis is by the linearized slender-body theory that has been 
widely used for steady flow phenomena. Although its reliability is 
limited to smooth slender bodies^ it is the most accessible approach to 
the problem that is already complicated by the necessity for considering 
time dependence. 
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The motion of the body is defined as the deviation from a figure 
of revolution originally parallel to the line of flight. Both rigid 
movement and deformation are represented by lateral translation of the 
circular cross sections expressed in terms of the corresponding 
translation of a deforming center line . 

The analysis is best suited for application to a pointed body. It 
can be formally applied to a blunt body but the approximate solution 
obtained is not usable near the nose. The type of shock wave formed by 
a blunt nose will alter the flow over the body and will affect the 
aerodynamic response to lateral motion. 

The solution is obtained in a definite integral that has not been 
found to be integrable in finite form. Evaluation of the integral has 
been accomplished by numerical methods. 

When the body radius is assumed to be vanishingly small, the solution 
degenerates to the well-known apparent mass approximation (ref. 2) that 
is obtained by regarding any cross section as a section from a two- 
dimensional incompressible flow. This approximation has the special 
feature that the aerodynamic force on any cross section depends only on 
the cross section and downwash, and their derivatives at that section. 


NOTATION 


a speed of sound 

b center of rotation for rigid body motion 

c diameter of base 

Cm pitching-moment coefficient, m , om . er ft , 

<3co S5 

I moment of inertia of equivalent air mass 

k reduced frequency, — 

a 

Z length of tody 

AL local lift 

ALp real part of 

ALj imaginary part of e^^AL 

M Mach number 
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<1 

Qoo 


angular velocity due to pitching 

free- stream dynamic pressure, i pV 2 

d 


(r,0 ,x) cylindrical space coordinates 
R radius of body 

Re real part 

S area of base 

t time 

v volume of body 

V free- stream velocity 

W downvash 

W transformed downvash 

W G amplitude of sinusoidal downvash 

Z. displacement of center line 

a angle of attack 

p air density 

cp velocity potential 

cp transformed velocity potential 

t -i- 

sin 6 

w circular frequency 

When a, a, and q are used as subscripts, a dimensionless derivative 
is indicated. 


r - 


Si + % " 


Sc, 


■m 


bO 


■m 


_a(ds/v) a(qc/V)J 
A dot above a symbol denotes a derivative with respect to time. 
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STATEMENT OF THE PROBLEM 

The linearized differential equation of the velocity potential cp 
for time dependent supersonic flow is, in cylindrical coordinates, 

-U 2 <Pxx + <Prr + i <Pr + ^ Use - % %ct - pr = 0 U) 

The coordinates are fixed in space and the undisturbed flow ahead of the 
body is of velocity V in the direction of the positive x axis. The 
derivation of this equation and a discussion of its limitations are given 
in chapter I of reference 1. The form given here is equation A 2 of 
table I of that chapter,, transformed to cylindrical coordinates. The 
axial distance downstream from the nose is given by x, and r and Q are 
polar coordinates in any plane of constant x (see sketch (a)). The 
radius of the cross section is R(x) and the displacement of the line 
through the centers of the cross sections is Z(x,t) in the direction 
0 * jt/2. 

Then in the linearized approxi- 
mation the boundary conditions for 
the part of cp that produces lift 
are 

I>r]r=R(x) = W(x,t) sin 9 (2) 

cp 0 as r ■* co (3) 

where 

W = Zt + VZx (4) 

The displacement of the center line Z(x,t) can represent any plane 
translation, rotation, or deformation of the body. The variable W(x,t) 
can be regarded as a form of downvash function if 9 = s t / 2 is directed 
downward. The solution satisfying boundary condition ( 2 ) in the slender- 
body approximation is dependent only on the function R 2 (x)w(x,t) for 
any W obtained by equation (4) , 

The local pressure change on the body is given by 

Ap » -p[cp t + Vcp x ] + . . . (5) 

and the local lift per unit length is 

p 

AL = / R Ap sin 9 d6 (6) 

Jo 

This lift AL is positive upward when Z and W are defined positive 
downward . 



Sketch (a) 
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Quadratic terms that are sometimes important in Ap do not contri- 
bute to AL. Also, the axially symmetric disturbance velocity component 
of the flow does not produce lift and need not be considered. 


ANALYSIS 


The solution of equation (l) with boundary conditions (2) and ( 3 ) 
is obtained in appendix A. The surface boundary condition is applied 
through the usual slender-body approximation yielding a distribution of 
sources on the center line. The source strength at any point depends 
only on the body radius and motion at that cross section. 

The method of solution is based on the use of Laplace transformations 
in t and x. The time transforms of cp and W are given in equations 
(Al) and (A2) of the appendix. The solution is obtained from the time 
transforms for the transform of AL in equation (Al8) and the Laplace 
inversion for AL Is given in equation (A19) . This is a general result 
for any distributions of body cross section and center~line motion 
allowable in the slender-body approximation. 

A simplified formulation is obtained if the time dependence is 
sinusoidal as defined in equation (A20) . Then a steady sinusoidal 
solution for AL is obtained by equations (A17) and (A2l) . This is the 
form used in the subsequent calculations for the results of the slender- 
body theory. 

A further simplification is obtained by the known apparent mass 
approximation that is formed here as a limit form for small body radius 
as stated in equation (A22) . This approximation is used for comparison 
to the results of the more exact slender-body theory. 


In the computations b only sinusoidal time dependence will be considered, 
that is W(x,t) = W 0 (x)e^ a) *^. 


We assume that W 0 (x)R 2 (x) is continuous over the length of the body. 
Since, in general, W 0 (x)R 2 (x) cannot be represented by the same expression 
over the entire length of the body (e.g., parabolic arc nose and cylindri- 
cal body) , we divide the body into k sections such that for each 


x j < x ^ x j+i 


W 0 (x)R 2 (x) = Wj(x)Rj 2 (x) 

Then AL is given "by equation (B7) in appendix B which is a deri- 
vation of the computation forms. Equation (B7) is sufficiently general 
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to cover a wide range of "bodies and downwashes. However, we note that 
unless the first derivative of W 0 (x)R 2 (x) is continuous at x=x^ , a 
singularity occurs at x=xj + (3r. This is where the Mach cone originating 
at x=Xj, r=0 intersects the surface of the body. Physically, there is 
no reason for such singularities. Mathematically, they exist because, 
while the problem prescribes boundary conditions on the surface of the 
body, the slender-body approximation gives a solution for an infinitely 
thin body, that is, the boundary conditions on the axis. Hence, linearizec 
slender-body theory as used here breaks down locally for these cases. 

The problem can be circumvented by approximating the corners by polynomial 
curves in such a manner that (d/dx)W 0 (x)R 2 (x) is continuous along the 
entire length of the body. This procedure is unnecessary unless some 
desired x falls very close to x=x j + J3r. Even with the above modifi- 
cation we do not know whether reasonable agreement with experiment would 
result for bodies with discontinuities. The best results from slender- 
body theory are expected for smooth bodies. 

The computation forms have been programed for the IBM 704- for cases 
where Wj(x)Rj 2 (x) can be expressed as a polynomial of degree < 9; a nd 
J < 7 • Thus the program can be used directly for certain deforming 
bodies and is not confined to rigid body motions. The program also 
includes a numerical integration of the local lift to obtain the 
generalized forces defined by 


F« = 


ALqXx) dx 


(T) 


where q(x) is a polynomial in x, and 2 is the total length of the 
body. The machine computing run time varies from four minutes for a 
simple case where J=1 and 20 values of x are calculated to 30 minutes 
for a more complicated body with J=7 and x calculated for 52 values. 


RIGID BODY MOTIONS 


Both linearized theory and apparent mass theory were applied in the 
calculations of the forces on a simple body in rigid motion. Rigid 
body motions were chosen for the comparison because they were the simplest 
to compute and because experimental data can be obtained to check the 
theories much more readily than for flexible models . Rigid body oscilla- 
tions should be sufficient to indicate theoretical accuracy. Moreover, 
the results are directly usable in the dynamic stability problem of the 
rigid missile. 

For the case of a body of revolution pivoting about a point b, see 
sketch (b) , we have 


W.(x)e iwt 


= co 0 e 


icot 


[ iu(x-b) + V] 


(8) 
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We divide AL into its real and imaginary 
parts: 



where 


S area of the hase 


c diameter of the hase 


If equation (B7) is used in equations (10) and (ll) > the integrations 
must he done numerically since ALp and ALj cannot he integrated analyt- 
ically* However,, the machine program permits the choice of as many 
values of x as desired for AL. 


The apparent mass formulas for C m and + C m may he written: 

cl oc» q 


„ 2v 2 fn ,_\ , w 2 I 
S5 ■ 5 + 


( 12 ) 


v volume of the hody 

I moment of inertia of an equivalent air mass referred to the center 
of rotation 


+ Cm = 0 (l-hf (ref. 3) 


( 13 ) 


Note that Cm* + Cm is independent of Mach numher and frequency, 
u* q 
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COMPARISON AND DISCUSSION 


The coefficients C™ and Cm. + C m were calculated for a hody with 

Iu a q 

parabolic arc nose and fineness ratio 10 (see sketch (b) ) . The equations 
for the radius are: 

R(x) in = (30-x) 0 < x < 15 in. 

R(x)i n = 2.5 15 < x < 50 in. 

Both the linearized theory and apparent mass formulas were used. Also, 
wind-tunnel tests were run on such a “body in the Ames l4- Foot Transonic 
and Unitary Plan Wind Tunnels. Figures 1 and 2 show the results for 
C m ^ and C m ^ + C m ^ plotted against the center of rotation for Mach 

numbers 1.1 and 1.2. At Mach number 1.1 and frequency 10 cps (k=5rt/72) , 
both theories agreed well with experiment. The agreement at MN1.2 was 
less satisfactory. At low Mach numbers and frequencies there seems to 
be little reason for using the slender-body theory solutions instead of 
the simpler apparent mass ones, at least in the case of rigid hody 
motions. Figures 3 and A show the effect of increasing Mach number on 
the theoretical values of C m ^ and C m ^ + C m . Both slender-body and 

experimental results show increase in coefficients with increasing Mach 
number, hut the apparent mass results show little change. Therefore, it 
appears that for the higher Mach numbers, slender-body theory solutions 
are superior to those derived from apparent mass. Figures 5 an( i 6 show 
the ratio of results of the slender-body theory to those for the apparent 
mass as Mach number increases. In figure 5 th e curves for b Jl = 0.4 

and 0.6 practically coincided with b Jl = 0.5; hence they are not 

included. Figures 7 and 8 show the effect in increasing frequencies on 
theoretical values of C m ^ and Qm^ + Cmq/ The effect of frequency is 

about the same for both theories. The difference in the absolute values 
is primarily a Mach number effect. Figures 9 and 10 compare the local 

lift as computed by slender-body theory and apparent mass theory at 

M=1.2 and M=2, A significant increase in the difference in the lift 
distribution at the higher Mach number is noted. 


CONCLUDING REMARKS 


For low supersonic Mach number and low frequency oscillation, the 
linearized slender-body theory used here has no advantage over the 
conventional apparent mass approximation using the mass of displaced air. 
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For higher Mach numbers the theory shows the local lift and, in 
general, the integrated lift and moment, to he larger than those obtained 
by the apparent mass approximation. This increase has also been found 
experimentally . 

For higher frequency of oscillation the theory indicates an increase 
over the apparent mass result but this does not appear likely to be large 
for frequencies that might be encountered. No experimental results have 
been obtained for the effect of frequency. 

The slender-body theory is limited to smooth pointed bodies. Rapid 
but continuous changes of cross section even with continuous slope cause 
■uncertainties in the results. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Feb. 28, 1961 
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APPENDIX A 


SOLUTION OF THE LINEARIZED SIENDER-BODY EQUATIONS WITH TIME DEPENDENCE 


The differential equation (l) 


-fnw + 9^ + - <P„ + rs <P, 


2 V 


cp. . = 0 


xx T r rr T r Y r r ^ ^99 Y xt W y tt 

is to he solved for the boundary conditions (2) and (3) 

[Cp r ] i-R(x) = W ( x ^) sln 6 
cp -*• 0 as r -* oo 

The use of Laplace transforms provides a convenient method of 
solution. Time transforms of cp and W are taken in the form 


( 1 ) 


( 2 ) 

(3) 


cp — / e~ s ^cp dt 

o 


(Al) 


W = / e-stw dt, Re(s) > 0 

Jq 


(A2) 


The remaining variables are unaffected by this transformation and equation 
(l) becomes 


py jg — gS _ 

^xx + 'hr + 7 L + P? - S? \ - 55 * - 0 


— 9 . + -% 9 

r r r 2 D 

The boundary conditions ( 2 ) and (3) become 


[ ^ R (x) " W sin 0 


(A3) 


(AL) 


9 -*• 0 as r -* 00 (A5) 

The transformations take this form -when 9 and W vanish for t < 0 . 

The only solution of equation (l) compatible vith boundary conditions 
(2) is of the form 


cp(x,r,0,t) = \Jr(x,r,t) sin 0 


(A6) 


n 


and after transformation (Al) 

cp(x,r,e) = ^(x^r)sin d 


(A7) 


The differential equation and boundary conditions now become 

-B 2 \f +t !.{(?_ 2Vs £ - £ f = 0 

H v xx T rr r T r r 2 ^ a2 ’x a.2. * 


(AS) 


and 


^r^reR( x ) 


(A9) 


f -*• 0 as r -* oo (A10) 

Separation of variables in equation (A8) yields a solution in the 

form 


e^KiCtrr) (All) 

where Ki is the first order modified Bessel function of the second 
kind and 


with 


2l2 




2VsA 


a c 


Re(cr) > 0 for Re (A) > 0 and Re(s) > 0 


The selection of the Ki solution of the modified Bessel equation is 
necessary to satisfy boundary condition (A10) . 

The solution of equation (A9) is sought as an integral of solutions 
in the form of an inverse Laplace transform 


t(x,r) 


X: 


C+i» 


ioo 


^KiCor) f (A)dA , 


c > 0 


(A12) 


At this point the slender-body approximation is introduced by replacing 
Ki by its approximate value for small r 


Ki(cr) 


_1_ 

ar 


(A13) 
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and; to this approximation 


t(x,r) 


1 

2 jti 


C+loo 


J\X 


C~loo 


— f(A)dx 
err 


Differentiating and substituting in equation (A 9 ) 


W(x) = 




2jcL 


c+ioo 


Ax 1 


C-loo 


crR^x) 


f(A)dA 


OTj rearranging 


C-ico 


and hy the Laplace in-version formula 


= _ f e“^ x, R 2 (x')W(x’)dx' 


Putting this in equation (A12) 


p G+loo p 1 - 

\jr(x,r) - / e^KxCcJr)^ dA / R 2 (x ! )W(x ! ) dx f (Al4) 

2 rtiJ J 

C-loo o 


This is a repeated integral that does not permit reversal of the 
order of integration. It can he put In a form that does permit the 
reversal hy expression as a derivative. 


?(*>r) = ~ — I 

Sr 2jri J c _ ic 


C+Ico 


s^KoCdr) dX J' e~^ <x R 2 (x’ )w(x f )dx' 


a 

Sr 


00 

R 2 (x')W(x')ax’ 


1 

2jtL 


i c+it 


4 . 


s A(x-x' )K 0 (sr)dX 


(A15) 


ioo 


The integral in the brackets in equation (A15) can he reduced to a 
standard form hy the substitution 


Vs 

0 2 a 2 


5 = A + 


and it becomes 


e 


Vs 

j3 2 a 2 


(x-x' ) 


1 

2jfi 


r Cl+1 ^( x -x 

Ci-ioo 




with 


Cx > 


s 

P 2 a 


The integral is obtainable from (Vf) , page 28 A, reference 4, and is 
reducible to elementary form by noting that 


I-x /2 (*) 


cosh x 


Then, observing the reduced limit on the x' integration in (A15) 
Ms 

■ {x-x' ) r 

s 


t(x,r) = J 


x-Pr 


2 (x-x') 

e H cosh 


.p s a 


JJx-x' ) 2 -P a r 2 


'/(x-x' ) 2 -p 2 r 2 


R 2 (x' )W(x' ) dx' 
(Al6) 


From equations (2) and ( 5 ) 


r2K 

AL = -p / R[cp-t + V<p x ]sin 0 d0 


J o 


and from equations (6) and (A6) 


nsit 


AL = ~p / R[t t + Vitf x ]sin 2 0 d6 


J o 


= -*pR[V t + W X 1 

Then the time transform of AL is 


(AIT) 
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For the general case, with Z = constant and therefore W = 0 for 
t < 0, the transform W is obtained hy equation (A2) and then f is 
formed by the Integration in equation (Al6) . Then after substitution 
in equation (Al8) an inverse Laplace transformation yields 


AT i r c+io ° st - . 

^ = ’2rti J e AL ds , 

c— ioo 


c > 0 


(A19) 


For a sinusoidal motion of the form described by 


Z(x,t) = Z n (x) 


iwt 


(A20) 


W(x,t) = W 0 (x)e lut 
it is possible to obtain a steady oscillatory solution. In this case 

W 0 (x) 


¥(x) = 


s-iw 


and using equation (Al6) in the Laplace inversion 

^ C+lOQ 


f 


1 

2jri 


e st f ds 


C-loo 


and reversing the order of Integration 


_ r x " Pr r 2 (x’)w 0 (x') _i_ p 

O J~f , \ 2 a 2 ' 2 2?ti -4 

v (x-x ! ) -P r 


C+ico si 

e 

■c-ioo 


M / % 

* - ^ (x_x,) 


cosh 


s 

_P 2 a 


/(x-x') a -[3 a r 2 


ds 

s-iu 


For t > x/V-a the oscillation of \|r is steady and the inner integral 
is given by the residue at s=iu and 


t = e iu)t A 

<5l\ 




io>M / f s 

v x ~ x ) 


cos 


0) 

J3 2 a 


— — 
7(x-x') -p r 


4 


,%2 o 2 2 

x-x' ) -p r 


R^x’ )W D (x' ) dx' 


(A21) 


From this result AL is formed by substitution in equation (A17) . 
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The above results were obtained by the usual slender-body approxima- 
tion. The same approximation can be used to obtain a simplified result 
from equation (Al4) . If the approximation (AI 3 ) is introduced, equation 
(Al4) becomes 


e^dAj^ e " XX ' R2 ( x '^^ x, ^ dx 

C-ioo 


= - i R 2 (x)W(x) for x > 0 
r 


from the Laplace transformation and its immediate inversion. A time 
inversion of and $ yields 

i|r =s - — R 2 (x) W(x ,t) 
r 


and 



R 2 (x)w(x,t) 


(A22) 


This is the usual time dependent generalization of the apparent mass 
concept as used by Munk (ref. 2) . It should be noted that AL of 
equation (A22) depends on R and ¥ and their derivatives at the point 
x at time t and is not affected by the shape of the remainder of the 
body or by any events at earlier time. 



l6 


APPENDIX B 

DERIVATION OF COMPUTATION FORMS FOR LINEARIZED SLENDER-BODY THEORY 


From equations (AI7) and (A2l) we have 

AL = jtpR[* t + Vt x 3 


t|t = e 


iut 3 


■>x-Pr 


if <“•> 


COB 


0) 

LP^a 


/(x-x- ) £ -p 2 r 2 


V (x-x' ) £ -0 2 r 2 


(A17) 

R a (x' )¥ 0 (x' ) dx' 
(A21) 


where and i|r x are evaluated at r=R(x) . 

Dividing the Body into k sections By planes perpendicular to the 
axis so that for each Xj < x < Xj +1 

W q (x)R 2 (x) = Wj(x)Rj 2 (x) 


we have 


♦ = 


j=0 X J 


iwM 


(x-x ! ) 


cos 


L 


a -(3 2 r 2 


Wj(x' )Rj 2 (x' ) dx' 


V(x-x’) 


2 „2 2 
-P r 


Xq = 0 

x J+l = x "Pr j 

(Bl) 


To carry the differentiation under the integral sign, we make the 
substitution: 


(x-x 1 ) = Pr sec 0 (B2) 

Applying (B2) and carrying out the indicated differentiation gives: 
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♦ = 


iwt 


— / \ 3^a 

(x-xj)e 


•z 

JaO 


iuiM , . 

Pi < x - x j> 


7(x-xj) 2 -p 2 r 2 


3 2 a 


^(x-x,) 2 -p 2 ^ 


w j (x j ) R j 2 ( x j ) 


x-X j+i) e 


it*)M t x 

pi ( x - x j+i) 


J-l 

♦P , 

>o sf(x-xj +1 ) 2 -3 2 r 2 

%) 


^•V(x-x J+1 ) 2 -P 2r2 J Wj(xj +1 )Rj 2 (x J+1 ) 


_ arc cos , 

J Vx 


-E'/ 


iwM 

FI 


3r sec 0 


J=o 


V X - X W 


IojM 

3a 


sec 2 0 cos fe 3r tan 0^) W^(x-3r sec 0 )Rj 2 (x- 


sec 0 tan 0 sin tan e )v x-3r sec 0)Rj 2 (x-3r sec 0) 


+ 3 sec 2 © cos (A* tan 0^) Wj(x-3r sec 0)Rj 2 (x-3r sec 0) J 


d0 


If ve assume that the W.(x)R^(x) are continuous at the x . y that 

J J J 

w j (x o+i )R j 2(x j+^ = W j+l (x O+i )R 0+l 2(x j+l ) 
and that R(o) = 0, the constant term vanishes. 


•3 0 arc cos 


arc cos te) . n* 


y - —sp 3r sec 0 


3 a 


\x-x J+1 y 


sec 0 cos ^5— 3r tan O^J Wj(x-3r sec 0 )Rj 2 (x- 


^ tan 6 sin (-£- 3r tan 0^) Wj(x-3r sec 0)Rj 2 (x-3r sec 0) 

■ 3 sec 0 cos (pi pr tan 0^ Wj(x-3r sec 0)Rj 2 (x-pr sec 0) J sec 0 d0 


From equation (A17) 


AL = -jrpR[tt + Vt x ] 


■3r sec 0) 

(B3) 

is, 

(■Bh) 

-3r sec 0) 

(B 5 ) 

(A 17 ) 


tt = 


(b6) 
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The calculation of is straightforward. If we apply (bU) to 

the result and change the integration variable by the substitution 

| = £$r tan 9 

Then, letting r « R(x) , we obtain the form 


AL = jtpe 


iut 


IuMk 

£ 2 a 


2Vx 


v/x 2 -^ 2 ! 2 


-7x a -P 2 r«j W o (0)H o (0)| 5 | r Rjx') 


iuiM 

J \X— X j} ^ ^ 

+ Z - — J ~ COS [w j (x')b/(x-)-W j . 1 ( x -)Hj. 1 2 (x-)1 

j=i 7 (x-x j) 2 -^ ^ J L Jx'-xj 


J J(x-Xj) 2 -p E r £ 
j=o 7(x-x J+1 ) B .p 2 r 2 


(- 5? 008 + 6ln ® V(x- 


+ (M^-t^cos 1^. + . . . Bin S w ( x _ 7| 2 +g 2 r 2 )E 2 ( x - J{ 2 -if) 2 r 2 




v cos ft S w J (x ' ^ )r j 2(x - ■/I^5 5 ?) 


(B7) 


For [ (d/dx)Ro(x) ] X _ Q < 00 (i.e., sharp-nosed bodies) , the first term 

vanishes. If the slopes are continuous at the x. as well as the down- 
wash and its first x derivative, all the constant terms are zero. 
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APPENDIX C 


DERIVATION OF AND C mi + C m ^ FROM APPARENT MASS THEORY 


From equation (A8) we have 

AL = jtp + V W(x,t)R 2 (x) 

In the case of sinusoidal time dependent motion, 

W(x,t) = W 0 (x)e iut 


(Cl) 


(C2) 


For rigid tody motions 


where z = a 0 (x-h) e'^*' and h is the center of rotation, 
in (Cl) 

AL = itpa 0 ^~ + V [ ico(x-b)+V]R 2 (x) e iu>t 


W(x,t) = (zt + Vz x ) 


(C3) 

Substituting 

(C4) 


AL = rtpe^ W " t a 0 -T-u 2 (x-'b)R 2 (x) + V 2 R 2 (x) + icoV 


d/_ 

dx 


(x-b)R 2 (x) + R 2 (x) 


(C5) 


ALr = a Jtp 


-to 2 (x-h)R 2 (x) + V 2 -iL R 2 (x) 

dx 


(C5a) 


ALj = jtpooVa c 


— (x-h)R 2 (x) + R 2 (x) 
dx 


(C5b) 


Substituting (C5a) Into equation (lO) and (Cpb) into equation (ll) 
c m a * | -0 p(x-‘b) 2 R 2 (x) + (h-x) A R2( x ) & 


(C6) 
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The first term equals (w2/q 0 Sc)l where I is the moment of 
inertia of an equl valent air mass with respect to the center of rotation. 
Integrating the second term by parts gives 

|| (b-Z)R 2 (l) +§| f R^xjdx 


which equals [-2(Z-b)/c ]+[2v/S£ ] where v is the volume of the body 
Therefore 


C m = — - - (l- b) + 
m a Sc 5 q S5 


(C7) 


The solution represented by (C?) was not found in any of the 
available references* 




(b-x) ^ (x-b)R 2 (x) + (b-x) R a (x) 


dx (C8) 


Integrating the first term by parts and combining gives 

• If 


§T 


(x~b)R 2 (x) + (b-x)R 2 (x) 


dx 


(C9) 


Therefore 


Cm. + C m - — 


"a 


(CIO) 


(Equation (CIO) is derived by Tobak, et al., in ref. 3* This derivation 
is included here for completeness .) 
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Figure 1.- Experimental and theoretical for a = 0°, k 

for slender body with parabolic arc nose cone. 
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Figure 2.- Experimental and theoretical Cm a + f° r slender tody 

with paratolic arc nose cone for a, - 0° k = 5 ^/ 72 * 
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Figure 3 »- Effect of increasing Mach number on C ni( 
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Figure 8.- Effect of increasing frequency on Cm^ + Cun* 




31 



0 .2 .4 .6 .8 1.0 

X /I 

Parabolic arc nose cone 
b//=.5, k = 5ir/72 


Figure 9*~ Comparison of real part of local lift for M = 1.2 and M = 2. 
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Figure 10.- Comparison of imaginary part of local lift for M = 

and M = 2.0. 
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